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Abstract 

The null geodesic equations that describe motion of photons in Kerr 
spacetime are solved exactly in the presence of the cosmological constant 
A. The exact solution for the deflection angle for generic light orbits (i.e. 
non-polar, non-equatorial) is calculated in terms of the generalized hyper- 
geometric functions of Appell and Lauricella. 

We then consider the more involved issue in which the black hole 
acts as a 'gravitational lens'. The constructed Kerr black hole gravita- 
tional lens geometry consists of an observer and a source located far away 
and placed at arbitrary inclination with respect to black hole's equatorial 
plane. The resulting lens equations are solved elegantly in terms of Appell- 
Lauricella hypergeometric functions and the Weierstrafi elliptic function. 
We then, systematically, apply our closed form solutions for calculating 
the image and source positions of generic photon orbits that solve the 
lens equations and reach an observer located at various values of the po- 
lar angle for various values of the Kerr parameter and the first integrals 
of motion. In this framework, the magnification factors for generic or- 
bits are calculated in closed analytic form for the first time. The exercise 
is repeated with the appropriate modifications for the case of non-zero 
cosmological constant. 

1 Introduction 

The issue of the bending of light (and the associated phenomenon of gravita- 
tional lensing) from the gravitational field of a celectial body (planet, star, black 
hole, galaxy) has been a very active and fruitful area of research for fundamental 
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physics. The Munich astronomer Johann Georg von Soldner in 1801 pQ using 
Newtonian mechanics and assuming a corpuscular theory for light derived a 
value for the deflection angle in the Sun's gravitational field A<fi s = 0".875. 
Later, Einstein using the equations of general relativity [2J, derived a value 
of ~ 1".75 which is consistent with the findings of Eddington's solar eclipse 
experiment and subsequent measurements. 

Despite the importance of the gravitational bending of light, in unravelling 
the nature of the gravitational field and its cosmological implications not many 
exact analytic results for the deflection angle of light orbits from the gravita- 
tional field of important astrophysical objects are known in the literature. 

Recently, progress has been achieved [5] in obtaining the closed form (strong- 
field) solution for the deflection angle of an equatorial light ray in the Kerr 
gravitational field (spinning black hole, rotating mass). Thus, going beyond the 
corresponding calculation for the static gravitational field of a Schwarzschild 
black hole [3], [I]. More specifically, the closed form solution for the gravitational 
bending of light for an equatorial photon orbit in Kerr spacetime was derived 
and expressed elegantly in terms of Lauricclla's hypergeometric function Fr>- 
It was then applied, to calculate the deflection angle for various values of the 
impact parameter and the spin of the galactic centre black hole Sgr A*. The 
results exhibited clearly, the strong dependence of the gravitational bending 
of light, on the spin of the black hole for small values of impact parameter 
(frame dragging effects) [B] . In addition, in [5] , the exact solution for (unstable) 
spherical bound polar and non-polar photonic orbits was derived. However, the 
closed form analytic solution for the important class of generic (i.e. non-polar 
and non-equatorial) unbound light orbits were left out of the discussion in [6]. 

One of the unsolved related important problems so far was the full analytic 
treatment of the Kerr and Kerr-de Sitter black holes as gravitational lenses. 
The closed form solution of this problem is imperative since the Kerr black hole 
acts as a very strong gravitational lense and we may probe general relativity, 
through the phenomenon of the bending of light induced by the space time 
curvature of a spinning black hole, at the strong gravitational field regime. A 
completely unexplored region of paramount importance for fundamental physics 
and cosmology. 

It is therefore the purpose of the present paper to calculate the exact so- 
lution for the deflection angle for a generic photon orbit in the asymptotically 
flat Kerr spacetime therefore generalizing the results in [6] and solve in closed 
analytic form the more involved problem of treating the rotating black hole as a 
gravitational lense. The constructed Kerr black hole gravitational lens geometry 
consists of an observer and a source located far away and placed at arbitrary 
inclination with respect to black hole's equatorial plane. 

More specifically, we solve for the first time in closed analytic form, the 
resulting lens equations in Kerr geometry, in terms of the Wcierstrafi elliptic 
function p(z), equation (|85p . and in terms of generalized hypergeometric func- 
tions of Appell-Lauricella equations ([89j), ([82]). ((60)) . (|99|) . 

In addition, we calculate for the first time exactly the resulting magnifica- 
tion factors for generic light orbits in terms of the hypergeometric functions of 
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Appell and Lauricella. Our closed form form solutions for the source and im- 
age potitions of the lens equations and the corresponding magnification factors 
represent an important progress step in the extraction of the phcnomcnological 
and astrophysical implications of spinning black holes. 

The resulting theory is of paramount importance for the galactic centre 
studies given the strong experimental evidence we have from observation of 
stellar orbits and flares, that the Sagittarius A* region, at the galactic centre of 
Milky- Way, harbours a supermassive rotating black hole with mass of 4 million 
solar masses [7J, [SJ. Given the fact that the GRAVITY experiment [5] and 
the proposed 30 metre telescope (TMT) [7J aim at an accurary of 10 //arcs the 
images calculated in this work and formed near the event horizon of the spinning 
black hole should be a subject of experimental scrutiny. 

Previous efforts on the issue of gravitational lensing from a Kerr black hole 
were concentrated on various approximations as well as numerical techniques 
using formal integrals [TU ] . [TT ] . 

The material of this work is organized as follows: in sect ion EJ we present the 
null geodesies in a Kerr spacctime with a cosmological constant. In section [3] we 
describe the Kerr-lens geometry and relate the first integrals of motion to the 
observer's image plane coordinates. In section 0] we derive a formal expression 
for the magnification using the Jacobian that relates observer's image plane co- 
ordinates to the source position. This expression involves derivatives of the lens 
equations in Kerr geometry and in our contribution we shall calculate in closed 
form these derivatives in terms of the generalized hypergeometric functions of 
Appell-Lauricella. In section we derive constraints from the condition that a 
photon escapes to infinity and it is not caught in a (unstable) spherical orbit. 
These constraints on the Carter's constant and impact factor define a region 
usually called the shadow of the rotating black hole. For values of the initial 
conditions inside the region enclosed by the boundary of the shadow and the line 
with null value for Carter's constant there is no lensing effect since the photons 
cannot escape and reach an observer. We also discuss constraints arising from 
the polar motion. In section [7J we derive for the first time the closed form so- 
lution for the angular integrals involved in the gravitational Kerr lens, in terms 
of the generalized hypergeometric functions of Appell-Lauricella. The full exact 
solution for a light ray which originates from source's polar position and involves 
m— polar inversions before reaching the polar coordinate of the observer is de- 
rived: equations (|64| . (|60 |) . In section [Jj we perform the analytic computation 
of the radial integrals involved in the Kerr-lens in terms of the hypergeometric 
functions of Appell and Lauricella. In the same section, we derive the closed 
form solution for the source polar position in terms of the Wcicrstrafi elliptic 
function p(z, g^ g^) that implements the constraint that arises from the first lens 
equation ([?]). In sections 151 18. 3[ we apply our exact solutions for the calculation 
of the source and image positions for various values of the spin of the black 
hole and the first integrals of motion, for an equatorial observer and an observer 
located at a polar angle of 7r/3 respectively. We exhibit the image positions on 
the observer's image plane. In Appendix [S] we collect the definition and the in- 
tegral representation of Lauricella's multivariable hypergeometric function Fd ■ 
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In addition in appendix [XJ we prove in the form of Propositions, some math- 
ematical results concerning the transformation properties of the function Fjj 
which are used in the main text. 



2 Null geodesies in a Kerr-(anti) de Sitter black 
hole. 

Taking into account the contribution from the cosmological constant A, the 
generalization of the Kerr solution is described by the Kerr-de Sitter metric 
clement which in Boyer-Lindquist (BL) coordinates is given by |12j-[13|: 



d.s 2 



„ r „ (cdt — a sin 2 



2 P \ dr 2 P 2 ia2 



-d(T 



Ag sin 
S 2 p 2 
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We denote by a the rotation (Kerr) parameter and M denotes the mass of 
the spinning black hole. 

The relevant null geodesic differential equations for the calculation of the 
gravitational lensing effects (lens-equation) and for the calculation of the deflec- 
tion angle are: 



A(t> = 
where 
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We also derive the equation related to time-delay: 
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The parameters <E>, Q are associated to the first integrals of motion. The 
former is the impact parameter and the latter is related to the hidden first 
integral (due to the separation of variables in the corresponding Hamilton- Jacobi 
partial differential equation (PDE)). 



3 The Kerr black hole as a gravitational lens. 

3.1 Observer's image plane 

Assume without loss of generality that the observer's position is at (ro,0o,0)- 
Likewise, for the source we have (rs,0s,4>s) ■ We also assume in this section 
that A = 0. In the observer's reference frame, an incoming light ray is described 
by a parametric curve x(r), y(r), z(r), where r 2 = x 2 +y 2 + z 2 . For large r this the 
usual radial BL coordinate. At the location of the observer, the tangent vector to 
the parametric curve is given by: {dx/dr)\ ro x+(dy/dr)| ro y+(dz/dr)| r . z. This 
vector describes a straight line which intersects the (o;,/3) plane or observer's 
image plane as it is usually called [9]-|llj at see figffl 

The point {ai,(3j) is the point (— /3j cos 0o, cti, cos 6o) in the (x,y,z) 
system. Our purpose now is to relate the cti,/^ variables to the first integrals 
of motion Q. For this we need to use the equation of straight line in space. 
A straight line can be defined from a point P\(x\, y±, z\) on it and a vector 
£( e i,£2,£3) parallel to it. The analytic equations of straight line are then: 

x-xi _ y-yi _ z- zi ^ 
ei £2 £3 

Applying ^ we derive the equations: 



cos 9o — ro sin 6o 



cos t) — r cos t) 



r Q cos 6> §£ | r=ro + sin 9 Q r Q sm9 ^£\r=r cos(9 - r Q sin6» ^| r= 



Solving for o^, P i we obtain the equations: 

d6, 



&i = 



, dr \r=ra 

Now we have from the null geodesies that: 

e(M 1/2 



-I 
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Figure 1 : The Kerr black hole gravitational lens geometry. The reference frame 
is chosen so that, as seen from infinity, the black hole is rotating around the 
z-axis. 
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Using cqns (fT3)) , (fT^]) and assuming large observer's distance ro(i-e. ro — > 
oo) we derive simplified expressions relating the coordinates (a.;,/?,) on the ob- 
server's image plane to the integrals of motion 

$ ~ —an sin 9 (15) 
Q s /3 2 + (a 2 -a 2 )cos 2 (# ) (16) 

We can also express the position of the source on the observer's sky in terms 
of its coordinates (rg, 9s, 4>s) an d the observer coordinates. Indeed, the equation 
for a straight line can be determined by two points P\{x\, j/i, z\), P%{x2, V2, ^2): 

x-xi y-yi Z-Z-y ^ 



X 2 ~X\ J/2 - Vl z 2 - Z X 

Thus applying the above formula for the straight line connecting the observer 
and the source yields the equations: 

r r s sin6> s sin cj> s 

Ois = 

to — rs(cos9s cos Oo + sin#o sin0scos</> s ) 

-r r s (sin 9q cos9 s - sin 9 S cos (j) s cos9p) 

5 ro — rs (cos 9s cos 9o + sin 9q sin 9s cos </> s ) 

4 Magnification factors and positions of images. 

In the following sections, we shall perform a detailed novel calculation of the 
lens effect for the deflection of light produced by the gravitational field of a 
rotating (Kerr) black hole and a cosmological Kerr black hole (i.e. for non-zero 
cosmological constant A). 

The flux of an image of an infinitesimal source is the product of its surface 
brightness and the solid angle Aw it subtends on the sky. Since the former 
quantity is unchanged during light deflection, the ratio of the flux of a sufficiently 
small image to that of its corresponding source in the absence of the lens, is given 

by 

M= (A^ = m (19) 

where O-subscripts denote undcflected quantities [5] and J is the Jacobian of the 
transformation [x Sy Vs) ~> (a»,2/i)0- Writting x s = x s {xi,yi),ys = Vs{xi,yi) 
we can find expressions for the partial derivatives appearing in the Jacobian 
by differentiating equations ([J} and Indeed, the Jacobian is given by the 
expression: 

J = xw — zy (20) 



1 Recall in the small angles approxiamation: cti rg Xj, /jj « tq Vi- Also we define: 



x$ := — , VS ■= — 
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where we defined: x := ^r,y '■= ^§f-, z '■= 7^i w := ^r- Writting equations 
((U and flU) as follows: 



Ri{xi,yi) - Ai(xi,j/i,a;s,2/s,m) 
A0(xs,ys,n) - R 2 (x i ,y l ) - A2(xi,yi,x s ,ys,m) 

we set up the following system of equations: 



where 
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dA x 
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dA 2 
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dRi dA 1 

'l — ~ HI - ' P2 



dxi 



dxi 



dRi _ dM 

dyi dy l 



dR 2 8A 2 



dxi 



dxi 



dR 2 dAn 



dyi dyi 



Solving for x, y, z, w we obtain: 
/' 



1 



a,\a^ — a 2 a-i 



(26) 



The parameters n = 0,1,2,... and m = 0,1,2,... are the number of wind- 
ings around the z axis and the number of turning points in the polar coordinate 
8 respectively. We shall discuss the latter in detail in the section that follows. 



5 The boundary of the shadow of the rotating 
black hole and constraints on the parameter 
space. 

The condition for a photon to escape to infinity , which is also the condition for 
the spherical photon orbits in Kerr spacetime [B] , is given by the vanishing of the 
quartic polynomial R(r) and its first derivative (also in this case ^P| r=r , > 
0). Implementing these two conditions, expressions for the parameter $ and 
Carter's constant Q are obtained [6], [T4] : 



a 2GM i „2 m qGM„2 i „.3 r a ( -4a 2 + r ( ~ 3 9i M + r) 



, 2 



(27) 

The perturbed, from the radius r = r, Ilst , of unstable spherical null orbits 
in Kerr spacetime, and thus escaped photon, will be detected on the observer's 
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image plane, at the coordinates 



ro sm0 o a Jr^W) 



±v /_ r 3[ r ( r _ 3GM)2 _ 4a 2^M] _ 2a 2 r(2a 2 ^ + - 3r^^)z - a 4 (r - ^) 2 4 



ro sin o a(r-^) 

Equations (|28p were derived by plugging into equations (fTS)) . ([IB]) the values of 
the parameters Q, $ that correspond to the conditions for the photon to escape 
to infinity, equations (|27j) . A photon will be detected when the argument of the 
square root in eqn.([28|) is positive. In Eqn ([28|) . zq ■= cos 2 9o- 

With the aid of equations ([15]) and ([H>|) we derive: 

a 2 + f3 2 = <f> 2 + Q + a 2 z (29) 



Apart from the constraints expressed by equations (|27|) we also derive con- 
straints for the motion of light from the allowed polar region:^;,! < 0s,9o < 
^max- Indeed using the variable Zj := cos 2 6j, we have z m > zo where z m is the 
positive root of: 

-a 2 z 2 m + {a 2 -Q- $ 2 )z m + Q = 

Let us see how this can be understood. Defining: z m :— zo ~ x we derive 
the quadratic equation for x 

- a 2 x 2 - x(a 2 - Q - $ 2 - 2a 2 z ) - a 2 z 2 + z (a 2 - Q - $ 2 ) + Q = (30) 

with roots: 



£1,2 = 



a 2 + g + 2a 2 z + $ 2 T ^4a 2 Q + (-a 2 + g + $ 2 ) 2 
2a 2 



(a 2 + 0) - a 2 w + ^((a? + 0) - a^wof + 4a 2 0w o 



2a 2 



(31) 



where w<j '■= sin 2 9o- The "radius" <I> 2 + Q must be greater or equal than the 
boundary of the photon region defined by Eqs. (j27|) and the line Q = 0. The 
minimum of this value is reached when Q = and a — > l.The actual minimum 
value is ($ 2 (r) + Q(r)) min = 4. Thus, by Eq.® we have that a 2 + > 4, 
and since < a 2 wo < 1, it follows the inequality a 2 wo — {a 2 +01) < and 
consequently, x < 0. Thus we conclude that z m > zo- Similar arguments ensure 
that when zs > zo it follows z m > z$ |llj . 



(28) 



9 



6 Closed form solution for the angular integrals. 



Let us perform now the exact computation of the angular integrals which occur 
in the generic photon orbits in Kerr spacetime thereby generalizing the results 
of [BJ. In the case under investigation, we have to take into account the turn- 
ing points in the polar coordinate. A generic angular polar integral can be 
written: 

r02 /•max(zi,Z2) /■mm(zi,Z2) 

± / = / +[1 - sign(^ o e 2 )] / (32) 

J Q\ J min(zi .Z2) «/0 

where: 

61 o6 2 := cos6»iCos6» 2 (33) 



Indeed, using the variable 



z := cos 



we derive: 



1 dz 1 • as a a 

27 = z7T=z =slsn{ 2- d)dd 



(34) 



This is the result of the fact that in the interval < 9 < |, cos6> > and 
sin 8 > 0, while in the interval 5 < 6 < 7r, sind > 0, cos 9 < 0. The angular 
integration in the polar variable includes the terms: 



= ± 



± 



9 s 



± 



.in/, 



±...± 



On 



(35) 



The roots z m , z% (of Q(8) — ) are expressed in terms of the integrals of 



motion and the cosmological constant by the expressions: 



Q + $ 2 S 2 -H 2 ± y/(Q + $ 2 S 2 _ #2)2 + 4H 2Q 

z ™,3 = (36) 



and 



H 2 : =^[Q + ($_ a )2 S 2 ]+a 2 s 2 



For A = 0, the turning points take the form: 



a 2 - Q - $ 2 + ^4a 2 Q + {-a 2 + Q + $ 2 ) 2 



2a 2 



(37) 



(38) 



where the subscript "m" stands for "min/max". The corresponding angles are: 

(39) 



^min/max — ArCCOs(±-\/ Z m ) 



Now for 8j and 8 min i max in the same hemisphere: 



iln/i 



<W 



dz 



±Z/W) 2 M4 t/z(z m -z)(z-z 3 ) 



= h 



(40) 
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Let us now calculate the elliptic integral in eqn. (|40[) in closed analytic form. 
Applying the transformation: 



Z Z m ~\- £ ( Zj ^m) 



(41) 



our integral is calculated in closed form in terms of Appcll's generalized hypcr- 
gcometric function F\ of two variables: 



h = 



1 



) (\ 1 1 3 z m -Zj z m - Zj\ r(i)r(i) 



^1 77) 77) 77) 77) 



2|«l {/z m (z m - zj) V2' 2' 2' 2' z m 'z m -z 3 y r(3/2) 
On the other hand using the transformation: 



(42) 



(43) 



we calculate in closed form: 

1 f*i dz 



2 \ a \ Jo y/z(z m - z)(z - z 3 ) 
|a| z„ 



. 1 1 3 2j 2j(^m — 23) 



^1 1) 77) 77) 77) ) 



2 2 2 z m z rn (zj Z3) 



2 /^j (-Stti — 23) 

1 M z m (zj-z 3 ) (1113 z Tra z 3 (z m -z 3 ) z 3 (z m - 2:3) 



|a| $~z~ n 



Z3 



2 ' 2 ' 2 ' 2 ' z m - z 3 z m (zj - z 3 ) ' z m (zj - z 3 ) 



(44) 

In going from the second line to the third of (|4~4"|) we made use of the following 
identity of Appcll's first generalised hypergeometric function of two variables: 

Fi(a, P, p', 7, x, y) = (l- - yy-^F^ - a, (3, 7 - /3 - /?', 7, ^7, 2/) 

a; — 1 
(45) 

Likewise we derive the closed form solution for the following integral: 



1 

2M./o (l-z)t/z(z m -z)(z-z 3 ) 
Zj 1 Zj z m 1 



Z m |o| 1 - Zj fyZjiZj -z m )(z 3 -Zj) 

Fd L l _1 1 3 z 3 (l - z m ) Zj(z m - gg) 

\ 2 2 2 Z m {\ Zj) Zm Z m .{zj Z3) 



1 Zi 



-z 3 z 



1 



1 1 3 



1 2' 1 '2'2'2' Zj U 



(46) 



Producing the last line of equation (|46[) we used the following formula for the 
Lauricella function Fd '■ 
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Proposition 1 

F D (a,P,P',P", 7 ,x,y,z) = (1 - ^"^'(1 - x)^(l - z)~ p " x 



Fd [i — — P — p — P iP ii-, °r,y, 



1 > f i 

IE — 1 Z — 1 

Proof. Applying the transformation: 

(47) 

onto the integral: 
l 

7 — a — 1 / 



IR Fd = w a ^(l - «) 7 (1 - ^^(1 - 1/ y)^ (1 - dti (48) 

Jo 

we derive: 



1 — vy ) \ 1 — vy 



and thus we obtain the result: 



IR Fd = (1 - y r- a (l - x)-"{l - y)-?\l - z)-P" x 

x — 1 z — 1 

(50) 







F D (a,p,p',p",j,x,y,z) = (1 - y)^ Q ^'(l - x)^(l - z)^" x 
■ 

Likewise applying (|4"Tj) : 

$ f 2 " 1 dz 



' 0" o" x y z v 



x — 1 z — 1 



2|o| 4 (1 - z ) ^z(z m -z)(z-z 3 ) 

^ 2 / ( 2 m Zj ) 1 2 ^ /1^113Zj Z m Z m Zj Z m Zj 

2|a| y z m f/(z m - z 3 ) (1 - Zm) \2' '2'2'2'l-z m ' z m 'z m -z 3 

(51) 
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Let us now compute exactly the term: ± J^ max / min ; m (|35[) : 



± I =2 I 

Anin / max — Z m and U m [ n O $ max 

Equation (f5T|) for zj = 0, becomes : 

$ 1 2 /l 1 1 3 -z m 



2|«| ^(z m - z 3 ) (1 - ^m) \2' '2'2'2'l-2 m ' ' z m 



-3 



$ 1 1 ^ zr Z^ 1 1 1 1 

"TT^l n ! 1) n' X ' 



|a| 




z 3 ) (1 - 




$ 


1 


Vr 


a 


|a| 


V (Zm ~~ 


^3)2 




$ 


1 


7T 


i 


|a| 


\J (z m — 


z 3 ) 2 1 - 


- z 3 



2 2 1 z m z ra z% 

1 ( 1 ^3 ) 



2" 



2m ^3 ^3 



Z Z Z m 23 \ Z Z Z m Z3 Z m — Z3 



(53) 



On the other hand the angular integrals of the form ± f e ™ ln 7 max in equation ([5]) 
are solved in closed analytic form as follows: 



,in/max _ $ J {z m -z s ) 1 2 ^ / 1 ^ 1 1 3 z s - z m z m - z s z m - z s 



2|o| V z m \/z m - z 3 (1 - z TO ) \2''2'2'2'l-z m ' z m ' z m - z 3 

$ ZS 2S - 2m 1 



-[1 - sign(0 s o 9 mS ) 



\a\ z m l-z s f/z s (z s - z m )(z 3 - z s ) 



P I11 113 z s (l - z m ) z s z s (z m - z 3 ) . 

Fd 1,1,--,-,-, — -j r, — , — -, (54 

2 2 2 z m (l — zs) z m z m {z s -z 3) 
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An equivalent expression for the above integral is: 

1 



± 



2M 



{Zrn - Zs) 



V ' z m — Z 3 (1 — z m) 



F I I 1 I I 1 Zs ~ Z '' 
I 2' ' 2' 2' 2' 



-[1- S ign(0 s o0 ms )]i{/^ 2 ' 
a V 



1 



2S - ^ 



23 



, 1 x 1 _1 3 z s (l ~ ^3) z s (z m ~ z 3 ) 
' 2 ' ' 2 ' 2 ' 2 ' z s - z 3 ' z m (z s - z 3 )' zs - z 3 



2M 



(z m - z s ) 



1 



[l-sign(6» s o6' m s)] 



^ — Z 3 (1 — Zffi) 



1113 Z S -^ r 

1 2' ' 2' 2' 2' l-z m 



1 



ZS - Z 3 



2:i 



- 2:s 



L l-^3 
1 



1113 Z S (1-Z 3 ) ZS (Zm " ^3) 
b D I 2' ' «' «' «' 



2.S' 



2 ' 2 ' 2 ' z s - z 3 ' z m (z s - z 3 ) ' z s - z 3 
zs 



1 1 1 3 z s (z m - z 3 ) 



1 - z 3 \ 2 ' 2 ' 2 ' 2 ' z m (z s - z 3 ) ' z s - z 3 



In going from equation (|54[) to equation (|55[) we made use of the functional 
equation of Lauricella's hypergeometric function Fo , Proposition 2 (|127[) . and 
Proposition 3 which are proved in Appendix [X] 

Now, for a light trajectory that encounters m turning points (to > 1) in the 
polar motion we hav^: 



where: 




(56) 



-[1 - sign(0 s o 6 mS )} 



+ +[1 - sign(0 o o mO )] 



+2(m - 1) 



(57) 



m o ■= Arccos(sign(yi)- v /i^) = Arccos(sign(/3 i )- v /^), 



(58) 



2 Recall the constraints of section [5] 
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Ui is the possible position of the image and: 



m o, m odd 

7T — d m o 1 m even 



Thus we have that : 
A 2 (xi,y i ,xs,ys,'m) = 2(m - 1) x 



$ 2 (Zm- ZS) 

2\a\ V z m 



M ^(z m - z 3 ) (1 " *m) 2 
2 



(59) 



1 fl, 1,1,1, , 

^2 2 1 — z m z m 



i 



V 7 z m — Z 3 (1 — Z m ) 



~S 



, 1 x 1 1 3 z s - z m z m - z s z, 

2 2 2 2 1 Zui Z m Zfn ^3 



-[1 - sign(6»s o e mS )} 



l°l«m ^Z 5 (^S-2m)(^3-^s) 



1 1 1 1 3 z s (l - z m ) z s z s (z m - zs) 
2' 2' 2' z m (l - 25)' z m ' z m (z s - z 3 ) 



2|o| y z m 



1 



•V^^m — ^3 (1 — Zm) 



, 1 j 1 13 Zp - ~ gm - Zp 

2 2 2 2 1 2 m Zjn Zjn — Z3 



[1 - sign(6> o mO )] 
-Pi) fl.l,- 



$ z z - z„ 



l°l z ™ 1 " z o Z]z {z - z m )[z 3 - z ) 

113 Zp(l - Z m ) Zp Zp(z m - Z 3 y 

2' 2 ' 2' z m (l — zq) ' 2 m ' z m (zo — 23) 



Wc now calculate in closed form the angular term A\(xi, yi,xs,ys, Tii) which 
appears in equations (|2"Tj) ,p|. 

Indeed, the angular integrals of the form ± J e 9 ™" l/max , in equation (QJ, are 
computed in closed-analytic form in terms of Appell's generalized hypergeomet- 
ric function of two variables as follows: 



± 



'6s 



Ad_ 1 V(z m -z s ) /1 in z m - z s z m -z s \ r(i)r(i) 

V® 2\a\ Zjz m {z m -z 3 ) 1 U'2'2'2' z m ' z m - z 3 J r(3/2) 



signet's o 

1113 
^1 t;* 77 j 77 > 77: 



Zs(z,n-Z3) 
2m(zS-23) 



|a| ^z m ~ z 3 

2s(2m-23) Z S (z m -Z 3 ) 



2 ' 2 ' 2 ' 2 ' z„ 



23 Zm(zS - Z 3 ) ' Z m (z S ~ Z 3 ) 



(61) 
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Also the integral P0|) is calculated for zj = in terms of ordinary Gaufi's 
hypergeometric function: 

Jo $J z(z m - z)[z - z 3 ) 

= 2J ^?[^^"f(1^^^) (63) 
2\a\ \z m {z m -z 3 ) \2 2 z m -z 3 ) 

Thus we obtain, 



Ai(xi,yi,x s ,ys,m) = 2{m-l)-^—J Zm 7rF(-,-,l, 



2|a| V z m (z m — z 3 ) \2'2' ' z m - z 3 



i ^(^m - z s ) (\ i i 3 z m -z s 4i - z s \ r(i)r(i) 

=-f 1 



2|«l ^ m (z m - z 3 ) i V2'2'2'2' z m ' ^ - J T(3/2) 

2/ «s(Zm-2 3 ) 

+ [l-sign(^og mS )]^ Vz '" (Zs - 23) x 



|a| ^z m - z 3 

1113 Zs(Zm ~ gg) Zg(z m - Z 3 ) 

2' 2' 2' 2' z m - 23 2 m (zs - 23)' z m (zs - z 3 ) 



1 y/(z m - z ) (1 1 1 3 z m -z z m ~ Zo\ r (|)r(l) 
^1 



2|a| Vz m (z m -z 3 ) V2'2'2'2' z m z m - z 3 / T(3/2) 
+ [1 - sign(6< o mO )] 



|a| v'Zm - ^3 



I---- Zm z o( z ™~ z 3) z o(z m ~Z 3 ) \ 

\2' 2' 2' 2' z m - z 3 z m (z Q - z 3 ) ' z m (z - z 3 ) } 



For m = i.e. for no turning points in the polar coordinate the exact 
solutions for the angular integrals in equation © become 
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Ai(xi,yi,xs,ys) 



and 

A 2 (x i ,y l ,x s ,ys) 



± 



Ou 

9s 



+ (1 - sign(0 s Q 96)) 

-(1 - sign(0 s o So)) 

1 1 1 3 z„ 



2 l a l \/Zm(Zm - Z 3 )' 



m 



zi 



2' 2' 2' 2' 



2 M y/z m {z m ~ z 3 ) 



+ [1 - sign(6» s o 9o). 



,1113 Z m - 

ll 2'2'2'2' z r 



Z2 Zm - Z 2 



Z3 



zi(z m -z 3 ) 
z m {z 1 -z 3 - l 



1113 

XF 1 ^^,n>n> 



M \/z m - z 3 

Z m Z\(z m -Z 3 ) Zi(z m -Z 3 ) 



w _. _ — ^ >: 



Z 3 z m( z l — Z 3 ) ' Z m (z\ — Z 3 ) 



(65) 



2W 



(z m - Zl) 



1 



V 2 ™ z 3 1 — z »i 
Zm Z m 



1 1 1 3 Zl 
xF D [ -, 1, - , 

2 2 2 2 1 ^7T7 



Z77 



-■-■■/J 



2R 



^2 



1 



. 1 1 1 3 z 2 - 2 
XFd ^ 2' 1 '2'2'2'T^ 



m %m ^2 ^ri 



Zi 



Zl 



+ [l-sign(0 s o# o )] 



Z 3 



" ^3 
1 



Zl — Z 3 y'Zm — Z 3 
Z3 Zl 



/ll 13 gi(l - z 3 ) zi z, 

77 ' 77 ' 77' 77 ' > > 

\ 2 2 2 2 zi - z 3 z m zi - z 3 zi - z 3 

(66) 

where zi := min(zs,zo), z 2 := max(zs,zo)- 

Equations . (jBT)]) for m> 1 turning points and (|^5|) . (f{J^l) for m = turning 
points constitute our exact results for the angular integrals which appear in 
(|2"Tj) for the case of vanishing cosmological constant A. It is time to turn our 
attention to the exact computation of the radial integrals which appear in the 
lens-equations of the Kerr black hole. 
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7 Closed form solution for the radial integrals. 



We now perform the radial integration assuming A = 0. 

For an observer and a source located far away from the black hole, the 
relevant radial integrals can take the form: 




For instance, in the calculation of the azimuthial coordinate ([5]) the following 
radial integral is involved: 




where A := r 2 + a 2 — 2GM \ . In order to calculate the contribution to the 
deflection angle from the radial term we need to integrate the above equation 
from the distance of closest approach (e.g., from the maximum positive root of 
the quartic) to infinity. We denote the roots of the quartic polynomial R (eqn 
(|6|) for A = 0) by a, j3, 7, 6 : a > (3 > 7 > 6. We manipulate first the terms: 



a (r 2 + a 2 ) 

A y/R 



dr 



adr 

7r. 



2GM, 



2GM 



adr r°° a^M^dr 



R 



Ay/R 



It is enough to proceed with the term [§ 



(69) 



00 a 2GM r _ a 2 $ 

A^R 



dr 



(70) 



Expressing the roots of A as r+ , r_ , which are the radii of the event horizon 
and the inner or Cauchy horizon respectively, and using partial fractions we 
derive the expression: 



A^R 



-dr 



A 9 ° 



zdr 



A go 



zdr 



+ 



(r-r+)VR J a (r-r_)^i? 
A 9 ° 

a (r - r+) f/Jr - a){r - f3)(r - j)(r - S) 
A go 



zdr 



a (r — r_) yj{r — a)(r — (3)(r - yj(r - 5) 



zdr 



(71) 



3 The radial term 2 ^= is cancelled from the angular term: — J ^=|- 
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where A 9 ± are given by the equations 

AT = ±— c - (72) 

r + — r_ 

For polar orbits $ = and the coefficients in ([T2"|) reduce to those calculated 
in 0. 

We organize all roots in ascending order of magnitude as followqj, 

a M > a v > at > a p (73) 

where a M = a^+i, a„ = a M +2, ot p = and a* = z = 1,2,3 and we have 

that a,,-! > a^,—2 > a M _3. By applying the transformation 



or equivalently 



where 



UJZ — 1 



(76) 

we can bring our radial integrals into the familiar integral representation of Lau- 
ricella's Fd and Appcll's hypcrgcomctric function F± of three and two variables 
respectively. Indeed, we derive 



-^^(a^+i - 




dc 


H+ 




- Z){1- K 2 +Z )yi- [I 2 






zdz 


#+ 




- z)(\~ K%z)y\- y?-z 






dz 


H~ 




-z)(i-K 2 _z)t/i-n 2 ; 


A 9 _°uj 2 (a, l+1 - 


a M+2) 


zdz 


H- 


^(1- 


-Z){\~K 2 _Z)Z]\-11 2 Z 



77) 



where the moduli k \ , /i 2 are 





- <V 











(78) 



We have the correspondence a^+i = a, 0^+2 = /?, ap-i = r+ = a (1 _2, a M -3 = 7, a M = <5. 
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Also 



and = r±. By denning a new variable z' := wz we can express the con- 

tribution A0^° ,to the deflection angle, from the above radial terms in terms of 
Lauricclla's hypcrgeometric function Fjj 



AC 



-2A a ° s /uj{a ll+ x - a^ +2 ) 



H- 



Fr 



I I i I 3 1 , 2 ^ 

2'2' '2'2'w' +' M , 



AfySjap.)-! - a^ +2 ) 
-2Ai°- v /£7j(a M+ i - q m+2 ) 



.3115 1/2 /2 



r(3/2)r(i) 



r(5/2) 



H- 



,1113 1- /2 

Fd( 2'2' 1 '2'2'^'"-' A< 



H- 



.3115 1 /2 , 2 

Fd ( 2'2' 1 '2'2'^'' t -' /i 



r(3/2)r(i) 

r(5/2) 



(80) 



where the variables of the function Fd are given in terms of the roots of the 
quartic and the radii of the event and Cauchy horizons by the expressions 



1 



a, j 



a 



'fi+2 



5-/3 



ot, 



a 



a — r± 

fi-1 
a — 7 



(81) 
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An equivalent expression is as follows 

-2A 9 ^ y /uj(a l _ l+ i - a^ +2 ) 



AC 



H' 



L .. 1 3 1 /2 , 2 



A 9 ° Vu{a»+-L ~ cVkO / 1/1113 1 



+ 



__ 72 F 1 o> n> n> o> 



2' 2' 2' 2' w 



-2A?°Vw(a M +i -^+a) / 1 1 13 1 /2 /2 
i?D U'2' 1 '2'2'^' K -' M 



1 



1 



11 13 1 



1113 1 



K '2 Fl { 2' 2' 2' 2' 



Ri{xi,yi) 



(82) 



In going from to (f8"2"|) we used the identity proven in [B], eqn.(52) in [5J. 
Finally, the term ^= ,is calculated in closed form in terms of Appell's 
first hypcrgcomctric function of two- variables : 



dr 
71 



(a 



"M-2) 



r(l/2) /l 1 1 3 a 1 
T(3/2) 1 \2' 2' 2' 2' x ' w 

1 ^fi, 4z2,4z4i(M) 



V(a-7)(a-5)r(3/2) V 2 ' 2 ' 2' 2' a - 7' <5 - a. 
We exploit further the lens-equations (|21[) . Indeed: 

1 / *Z /II \ /^s 



Jii(o; i ,» i )-2(m-l) 



Inverting: 



Zm 7 — i ( 1 

— -irr — ,— ,1, 

2\a\ V z m (z m - z 3 ) \2 2 z m - 23 



VS? - ff 2 £ " .93 
(84) 



£s = P 



while: 



2(m- 1) — 



-^('i,i,l, 



2|a| Y z m (z m - z 3 ) \2 2 ' z m - z 3 
<t> s = R2{xi,Vi) + A 2 (x i ,yi,xs i ys,m) 



+ e 
(85) 

(86) 
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where p(z) denotes the WeierstraB elliptic function (which is also a meromorphic 
Jacobi modular form of weight 2 ) and the WeierstraB invariants are given in 
terms of the initial conditions by: 

92 = ^{a + pf-Q^ (87) 

93 = —(a + f3) 3 - Q— - Q— (88) 
ys 216 v H; 48 48 v ' 

Also a := —a 2 , j5 := Q+$ 2 , z$ = — _ a /4 arj d e is a constant of intergration. 

To recapitulate our exact solutions of the lens equations (|21|) are given by: 

1 , a / \ 2 r(l/2)_ /l 1 1 3 

=dr = yi, x s , ys, m) ^ ^wo/n^ i 



VR iV ' N /( a _ 7 )( a _ ( 5)r(3/2) ii V2'2'2'2'a 



2|a| V z m (z m - z 3 ) V 2 ^' ' z m - z 3 



v 7 (z m - is) c / 1 i i 3 z m - z s z r „ - z s \ r(i)r(i) 



M Vz m {z m -z z ) V2' 2' 2' 2' z m 'z m -z 3 y T(3/2) 

2/ Zg(Zm-^3) 

+ [1 - sign(6» s o toS ) — x 



1113 -Zm -Zs(z m ~ Z 3 ) ^S(^m ~ ^3) 

2' 2' 2' 2' z m - z 3 z m (z s - z 3 ) ' z m (z s - z 3 ) 



1 \/(^m - z ) /1 1 1 3 z m - z z m - z \ r(|)r(i) 

2|a| ^z m (z m -z 3 ) 1 U'2'2'2' z m ' z m - z J T(3/2) 

2 / ^0(^771-^3) 

1 M la a M V 2 ™( z o-z3) 

+ [1 - sign(6» o 9 m0 )\ 



\a\ y/z m - Z3 



'1113 z m z {z m -z 3 ) z (z m -z 3 ) 



2' 2' 2' 2' z m - z 3 z m (z - z 3 ) ' z m (z - z 3 )_ 
-<j) s = R 2 (x l ,y l ) + A 2 (x i ,y l ,x s ,ys,m) (90) 

and equation (|55|) . 

In the subsequent sections we shall apply our exact solutions for the lens 
equations in the Kerr geometry expressed by (|89)) , (|85t . (|90)) to particular cases 
which include: a) an equatorial observer: 60 — 7r/2 =>■ zo = b) a generic 
observer located at Qq = tt/3 => Zo = j. 
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8 Positions of images, source and resulting mag- 
nifications for an equatorial observer in a Kerr 
black hole. 



In this case (6q = 7r/2) , equations (p~5jl . (|16|l . become: 



Q 



-O-i sin vq = —o>i 



p\ + (a? - a 2 ) cos 2 9 = p\ 



Thus the length of the vector on the observer's image plane equals to: 



* 2 + p* = V$ 2 + Q 

Furthermore, we derive the equations: 



as rs sin Us sin <p s 

ro r Q - r s sin S cos 4> s 



Vs ■= — = 



Ps 
ro 



—rs cos #s 



ro - rs sin 9s cos 4> s 



or equivalently: 



as 
Ps 



tan 9s sin <f> s 



(91) 
(92) 

(93) 
(94) 

(95) 



8.1 Solution of the lens equation and the computation of 

We now describe how we solve the lens equations (|2ip using the properties of 
the Weierstrafi Jacobi modular form p(z) equation f|85[) and the computation 
of the radial and angular integrals in terms of Appell-Lauricella hypergeometric 
functions equations (|53")) , (j52")l . ([M]) . respectively. 

For a choice of initial conditions a, $, Q we determine values for the observer 
image plane coordinates a%,P^ see equation (|9ip . Subsequently we determine 
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the value of zs and therefore of Os that satisfies the equation 
f°° 1 

J a VR via 



2 T(l/2) (I 1 1 3 0-7 5-/3' 



V(a- 7 )(a-(5)r(3/2) V 2 2 ' 2 ' 2 ' « " 7' S - a 



2\a\ y z m (z m - z 3 ) \2'2' z m - z 3 

i \J{z m - z s ) /i i i 3 z m -z s z m - z s \ r(i)r(i) 

= i< l 77 1 7' 7' 77' ' 



2 M l/z m (z m -z 3 ) 1 V 2 ' 2 ' 2 ' 2 ' z m 'z m -z 3 J T(3/2) 

2/ zs(Zm-Z3) 

+ 1 - sign(0 s o fl mg ) — JL- x 
|o| yz m - z 3 

F (-- -- Zm ZS ( Zm ~ Z ^ Z s( z m ~ Z 3 ) 

\2' 2' 2' 2' z m - z 3 z m (z5 - z 3 ) ' z m (>s - 2 3 ) 

1 n -Fdu «■ 



2 M </(z m - zT) V 2 ' 2 ' «m-*3 



(96) 



using our exact solution for zs in terms of the Weierstrafi elliptic function 
equation (J85J). For this we need to know at which regions of the fundamental 
period parallelogram the Weierstrafi function takes real and negative values. 
Indeed, the function of Weierstrafi takes the required values at the points: x = 
j + uj', I G K of the fundamental region (p( f + uj'; 32,53) G R~)- Thus as the 
parameter I varies we determine the value of zs that satisfies equations (|85[) . 
(|96|) . The quantities uj, uj' denote the Weierstrafi half-periods. In the case under 
investigation w is a real half-period while uj' is pure imaginary. For positive 
discriminant A c = g\ — 27 ' g\, all three roots ei, e2, e 3 of 4z 3 — 32Z — g 3 are real 
and if the e% are ordered so that e\ > e-i > e 3 we can choose the half-periods as 

dt ' ^ dt (97) 



1 V u3 - 9it - g 3 ' J-00 V~4£ 3 + ff 2 i + .93 

The period ratio t is defined by t = u/ /uj. An alternative expression for the real 
half-period uj of the Weierstrafi elliptic function is given by the hypergeometric 
function of Gaufi^|: 

1 -->(iU e ^) W 



V e l- e 3 2 V 2 ' 2 ' ' 61-63, 

Having determined 6s by the procedure we just described we determine the 



5 Which is the closed form solution of the radial and angular integrals of the first of lens 
equations in eqn 1211 . 

6 The three roots are given in terms of the first integrals of motion by the expressions: 
ei = ^(-a 2 + Q + * 2 +3 V / 4a 2 Q + (-a 2 + Q + * 2 ) 2 ), e 2 = i (a 2 - S - * 2 ), e 3 = ^j(-a 2 + 

Q + * 2 - 3 v /4a 2 Q+(-a 2 + Q + * 2 ) 2 ) 
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n — fi O — OA PiA^Pi'i <b — 9 71 Q1 1 
u, — u.u, ^ — z^.u^tuuo, — — z. ( lyiiu 


n — fl O — 1 98 <E> — 3 R^Q 


rv ( GM \ 
«« 1-73-J 

o (GM\ 

Pi K -W ) 


9 71 01 1 
z. / iy i iu 


— o.ooy 


— ^.yo^^ooozoy 


' : ! c i777n87fi' : !QQ 


( 2 GM \ 
1 \ r C J 


1.359555 


— 1.9195 




-2.48221826 


0.178885 


m 


3 


3 








zs 


0.3161007914992452 


0.0026145818604 


Os 


55.79° 


87.069° 


A0(rad) 


-11.086 


7.09441 


05 


95.1794° 


133.52° 




0.5545341990201503500 


0.824718843878947 


uf 


1.3278669366032567973i 


2.9400828459149726i 



Table 1: Solution of the lens equations in Kerr geometry and the predictions 
for the source and image positions for an observer at Qo = tt/2,<Pq = 0. The 
number of turning points in the polar variable is three. The values for the Kerr 
parameter and the impact factor 4> are in units of while those of Carter's 
constant Q are in units of G f . 



azimuthial position of the source </>g by the second equation of (|2 1 [) : 

- (j) s = R 2 (x l ,y l ) + A 2 (x i ,y l ,xs,ys,m) (99) 

Let us give an example. For the choice Q = 24.64563^1^, $ = -2.719110-^, a = 
0.6^ we determine z s = 0.3161007914992452, m = 3 and A(j> = -11.086, <j> s = 
95.1794°. Keeping fixed the value of the Kerr parameter we solved the lens equa- 
tions for different values of Carter's constant Q and impact factor <!>. We exhibit 
our results in table [ffl 

Let us at this point, present a solution with a higher value for the Kerr 
parameter in Table [2] 

The positions of the images of Tables [T] and [2] on the observer's image plane 
arc displayed in fig J2] and fig J3] respectively. In the same figures the boundary of 
the shadow of the spinning black hole is also displayed. 

8.2 Closed form calculation for the magnifications. 

Wc outline in this subsection, the closed-form calculation, of the resulting mag- 
nification factors. It turns out that the derivatives involved in the expression 
for the magnification are elegantly computed using the beautiful property of 

^Assuming that the galactic centre region, SgrA* , is a Kerr black hole with mass: Mbh = 
lO 6 M0and a distance from the observer to the galactic centre: tq = 8Kpc, the second 
solution in Table [TJvill require an angular resolution of 19.3102/iarcs . This is in the range of 
experimental accuracy for both the TMT and GRAVITY experiments. 
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Figure 2: The two images of Table [T] on the observer's image plane. The value 
of the Kerr parameter is a = 0.6-^p^, while the observer is located at 9<j = tt/2. 
With red is the image solution, first column of Table Q] and with green the image 
solution, second column of table [TJ 





a = 0.9939, Q = 27.0220588123, $ = -2.29885534 


ft (#) 

/ 2 gm\ 


2.29885534 
5.198274599547431 

1.14942767 


( ro c 2 ) 


2.5991372997737154 


m 


3 








0.01378435185109 


^5 


83.2575° 


A0(rad) 


-11.243 




104.177° 




0.5505433970950226 




1.1288708298860726 i 



Table 2: Solution of the lens equations in Kerr geometry and the predictions for 
the source and image positions for an observer at 8o = tt/2, (/)q = for a high 
value for the spin of the black hole. The number of turning points in the polar 
variable is three. The values for the Kerr parameter and the impact factor $ 
arc in units of while those of Carter's constant Q are in units of ¥ . 
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Figure 3: The lens solution of Table [5] as it will be detected on the observer 
image plane by an equatorial observer. The boundary of the shadow of the 
black hole is also exhibited. 



the hypergeometric functions: namely, that the derivatives of the hypergeo- 
metric functions of Appcll-Lauricclla are again hypergeometric functions of the 
same type with a different set of parameters. It is a powerful property of our 
formalism which we exploit to the full in what follows. 
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d^B dz s 
dz s dxs' 
$ 1 1 



2 1 a | z m (1 — z m ) ^Jz m — z 3 

F (- 1 1 I - Zs ~ Zm 

W '2' 2'2' 1-z™ 



2=S 



-1/2 



2* 2/ ( z m ~ Zff) 

2|a| V z m 



1 



F I 5 9 I I ^ 

1 2 '^2'2'2' 1-* 



^ z m 23 (1 — z m) 



,'3 3 15 z s - x 
1 2' '2'2'2' 1-z 



m 



zs z n 



Z 3 J 1 — Zn 



1 2' ' 2' 2' 2' 1-2 



^3 



l/'i- 



^3 



-1 



(1-Sigll(0so0 ms )(-1) 



1 ZS" ^r, 



1 



.z m 1 - z s y/ 'z s (z s - z m )(z 3 - zs) 
zs z 3 (z m - 3z s z m + 2z's) - z s (z m (2 - 4z 5 ) + + 3z s )) 



2(1 - z s ) 2 z s (z 3 - z s )y/z s {zs - z m )(z 3 - z s ) 
113 zs(l-z m ) z s z s (z m -z 3 ) 

r D J-i -L, — 7T5 7T> TTi 



2' 2' 2' z m (l - z s )' z m ' z m (z s - z 3) 
zs z s - z m 1 r 



z i« (1 ^ z s) V z s( z s - z m )(z 3 - zs) 

115 z s (l-z m ) z s z s {z m -z 3 ) 



F D 12, 2 
Fd (2, 1 



1 - Zr, 



2' 2' 2' z m (l - z s ) ' z m ' z m (z s - z 3 ) / z m (l - zs)' 
115 z s (l-z m ) z s z s {z m -z 3 )' 



' 2' 2' 2' z m (l - z s ) ' Zm z m (z s - z 3) J z m 

-13 5 z s (l-z m ) z s z s (z m - z 3 )\ / -z 3 (z m - z 3 ) 



F D [2,1 



2 ' 2' 2' z m (l - z s y Zm z m (z s -z 3) ) \ z m (z s - z 3 ) 



(100) 



r s sin flg cos s sin <ft g 

^ = (g00l)x -2cosgssin^ (r °" rsSinesCOS ^ )2 1 (101) 

OX S \ Jl 
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Now we calculate the term: 
we derive the expression: 



dzs 



. Indeed, calculating the derivatives w.r.t. z$ 



dz s 



Now: 



i r(i)r(i/2) / 

2\a\ r(3/2) ^ 

i r(i)r(i/2) , 

2\a\ T(3/2) V 
3 1 3 5 z T . 



2-\/ z m (z m — zsj-s/ z m — zs 



Fi 



1 1 1 3 Z r 



2' 2' 2' 2' 



[1 - sign(6> s o ms )} 



(z m - Z S ) 

— 1 — Z — T x 

i 

Zm %m ^3 

Zs{z m 



.3315 z m - 
1 2'2' 2'2'~ 



1 



2 2 2 2 ■ _/ f i 

i — zs z m — Z S 



Z3 



ZS Zn 

7 

-1 



ZS 
Z3 



%m %3 
Z3 



l 

2\a\ \z m (z s - z 3 ) 
zs(z m -z 3 ) z s (z 



(-z 3 )\/z m - z 3 
z m (zs ~ z 3 ) 2 

z 3 y 



2' 2' 2' 2' z m — zs z m (z s - z 3 ) ' z m (z s - z 3 ) 

3 3 15 z m z s {z m -z 3 ) z s {z m -z 3 ) 



Zs(z m -Z 3 ) 

z m (z s ~z 3 ) 



\a\ yjz m - Z 3 
/3 13 5 

^1 77' 77' 77> 77) 



i 7 ! 



2' 2' 2' 2' z m - z 3 z m {z s - 23)' z m (z s - z 3 ) ) \(z s - z 3 ) 2 

Zs(z m ~ Z 3 ) Zs{z m - Z 3 )\ ( {-Z 3 )(z m - Z 3 ) 



2 ' 2 ' 2 ' 2 z r , 



Z 3 Z m (z S - Z 3 ) ' Z m (zs - Z3) J V Zm(2S - 23) 



(102) 



Ql ^^ll^(fT02l) X |f£ = (fT(ax 



2 cos 63 sin 0s x 



r s sin &s cos ps sin q>s 
(rp-rs sin 6> s cos<ft s ) 2 

Jl 



and 



(103) 



a 2 = ^- = (HQ2J) x |^ = (H02| x [ -2 cos S sin 6> s x 
dys dy s \ 



(ro—rs sin 8s cos0 s ) 2 



•/l 



(104) 

While for the a 3 , 0,4 terms which contrbute to the expression for the magni- 
fication, equation (|26[) . we derive the expressions: 



a 3 



dA 2 



dx s dx s 



(rprs sin8s—r s cos<fi s ) 
(ro—rs sin 8s cos cf> s )' 2 

Jl 



r s sin us cos us sin <p s 



Cool) 



(105) 
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(106) 



where J\ denotes the Jacobian: 



Jl 



d(x s ,ys) 
d(6 s ,ct> s ) 



(107) 



and 



dxs 
d(j) s 



d8 s 
dx s 
d6 s 
dys 
d<t> s 



(r| sin 6> s cos 6 S sm(/) s )/((r - r s sm6 s cos(f> s ) 2 ) 
Ji 

-[rpr s sinflscos^g - r| sin 2 fl 5 ]/((r - r s sinflg cos(f> s ) 2 ) 

Ji 

(r r s sin 6 S - r| cos <f> s )/((r - r s sin6> 5 cos0 s ) 2 ) 
Ji 

r r s cos 6*s sin (f> s /((r ~ r s sin 9 S cos ^ s ) 2 ) 



(108) 



5ys 



Ji 



In producing the results exhibited in eqns (|100l) . (jl02p in our calculations 
for the magnification factors we made use of the important identity of Appell's 
hypergeometric function F\ : 



and its corresponding generalization for the fourth hypergeometric function of 
Lauricella. Similar calculations that we do not exhibit in this written account 
of our talk, lead to the derivation of the coefficients j3 1 , f3 2 , (3 3 , /3 4 in terms of 
the generalized hypergeometric functions of Appell-Lauricella. 

A phcnomenological analysis of our exact solutions for the magnifications in 
Kerr spacetime will be a subject of a separate publication [15] . 

8.3 Source and image positions for an observer located at 



In this case, the coordinates on the observers image plane are related to the first 
integrals of motions as follows: 



d m+n F 1 (a,p,p', 1 ,x,y) 



(a, m + n)((3, m)(/3', n) 
(7, m + n) 

F\ (a + m + n, /3 + m, (3 + n, 7 + m + n, x, y) 



(109) 



dx m dx n 





(110) 
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Furthermore, our solution for the first lens equation (|59"f takes the form: 



1 a At 2 £(1/2) /l 1 1 3 /?- 7 

sJR ^/( a _-y)( a _(5)r(3/2) V 2 2 2 2 a~7 <3-a. 



2(m - , Zm -^F ( 1 1 1, 



2\a\ Y z m (z m - z 3 ) \2 2 z m - z 3 

i ^(«m - *s) „ A i i 3 z m -z s z m -z s \ r(i)r(i) 



2 M ^z m (z m -z 3 ) 12'2'2'2' z m 'z m -z 3 y T(3/2) 

+ [1 - sign((9 s o mS ) — -L- x 
H V z m - z 3 

F (\ I I - Zm Z s( Z m ~ z s) Z S (z m - Z 3 ) \ 

1 \2' 2' 2' 2' z TO - z 3 z m (z s - z 3 ) ' z m (z s - z 3 ) J 

i ^m-i) /ni3«™-J z m - i \ r(i)r(i) 



2|«| ^Z m (z m -Z 3 ) V 2 ' 2 ' 2 ' 2 ' «m '^"V r(3/2) 



(l/4)(z m -s 3 ) 

M 1 V *m((l/4)-23) 



[1 - sign(- o mO )] 



3 \a\ Z/z m - z 3 



X 



p /l 1 1 3 z m (l/4)( ) (l/4)( 

1 \,2'2'2'2'z m -z 3 z m (l/4-z 3 ) ' z m ((l/4) - z 3 ) / ' 1 

Let us give now examples of solutions for the images and source positions of 
equations ([TTT j) . (|85j) . 

For the initial conditions Q =25.64563, a = 0.9939, $ = —3.11 we calculate 
the parameter z$ of the source latitude position to be: z$ = 0.09097820848 
(8s = 72.4447°) for m = 3. The position at the fundamental period parallelo- 
gram that provides the above value of zg as a solution of (|111[) , (|55|) for three 
turning points in the polar variable is located at: 1 299548O690017123 w '' wnere 
the fundamental half-periods of the Weierstrafi elliptic function p(x, 32, 53) were 
calculated to be: 

to = 0.52792338858688228, J = 1.119903617249492 i (112) 

Also the azimuthial position of the source was calculated to be using (J9H)) and 
the calculated value of z s : cf> s = 2.67231589rad = 153.112° = 551205'lS 



The first solution is shown on the image plane of the observer, FigJH We 
observe that the solution lies close to the boundary of the shadow of the black 
hole. 

8 A</> = B,2{xi,yi) + A2(xi,yi,x s ,y s ,m) was calculated to be: A.cj> = — 12.0971rad so that 
the photons perform more than one loop and a half around the black hole. 
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a = 0.9939, Q = 25.64563, $ = -3.11 


a = 0.52, Q = 23.64563, $ = -2.85 


OH (^) 
ft (#) 


3.591118674 


3.29089 


-4.7611506980 


-4.58320084657 


\ro c- J 


1.79556 


1.64545 




-2.38058 


-2.2916004 


m 


3 


3 








ZS 


0.09097820848 


0.5980171072414 


Os 


72.4447° 


39.3474° 


A0(rad) 


-12.0971 


-11.8577 


<f>S 


153.112° 


139.395° 




0.52792338858688228 


0.5571026427501503 


w' 


1.119903617249492 i 


1.389041935594241i 



Table 3: Solution of the lens equations in Kerr geometry and the predictions 
for the source and image positions for an observer at Qo = 7r/3, §q = 0. The 
number of turning points in the polar variable is three. The values for the Kerr 
parameter and the impact factor are in units of 



ft a=0.9939 




Si 



Figure 4: The solution, 1 st column of Table [3] as it will be detected on the 
observer image plane by an observer at 9o = tt/3, 4>q = 0- The boundary of the 
shadow of the black hole is also exhibited. 
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9 Exact solution of the angular integrals in the 
presence of the cosmological constant A. 

There has been a discussion in the literature as to whether or not the cosmolog- 
ical constant contributes to the gravitational lensing. However, the debate has 
been restricted to the Schwarzschild-de Sitter spacetime [TB], [17], [T8]. Let us 
discuss now the more general case of gravitational lensing in the Kerr-de Sitter 
spacetime. 

The generalized solution for the angular integral (|54|) in the presence of A is 
given by: 



± 



z m - z s 



hs 2 \ H \ C 1 - VZm) t/z m (z m - z s )(z m - z 3 ) 



$ „ / 1 . 1 1 1 3 rj(z s -z m ) z s -z m z m -z s z m -z s \Yy 2 



(1 - z m ) V 2 ' ' ' 2 ' 2 ' 2 ' 1- rjz m ' 1 - z m ' z m ' z m - z 3 J r (|) 
(I 1 1 3 i](zs - Z m ) Zm — Zs z m -2sU(i 



2' i '2'2'2' l-77z m ' z m ' Zm - z 3 >/ T (|) 
S 2 z s z m - z s 1 



(1 - sign(0 s o 6 mS )) 



\H\ z m 1 - T]z s N /z s (z s -z m )(z 3 -z s ) 



—aFn I 1,1, — , — , — , A, — — , it | H _Fn ( 1, 1, 1, — , — , — , A, v. — a 

' ' 2' 2' 2' ' z m J 1-zs I 2'2'2' ' ' z™ 



(113) 



where: 



a 2 A 



ZS Z r , 



Z3 



Zm Zs — Z 3 



Zm V 1 



(114) 



Also the integrals ± J e 



2 J Q contribute the term: 



2(m- 1) x 



S 2 $ 



2|ff| (1 - ryz m )(l - z m ) ^/4( Zm - z 3 ) 

x F fi 1 i I 1 2 ?/(-^m) ~ z m 1 
D \2' ' 'rm-riZm'l-Zm' ' 



z-i- 



1 



21^1 ^ 



*3, 



1 - ?7z l; 



«4 



^ 1 1 3 ?/(-z m ) 1 z m 
2 2 2 1 Tjz m z ra Z3 



(115) 



Notice that for A = this reduces to equation (|53|) . 



33 



9.1 Closed-form solution for the radial integrals in the 
presence of the cosmological constant A. 



Assume A > 0. We need to calculate radial integrals of the form: 



^((r 2 +a 2 )-a<i>))4L 

<-±r V A 



We use the technique of partial fractions from integral calculus: 
— 2 A 1 A 2 A 3 A 4 



A, 



-((r 2 +a 2 )-a$)) 



r — r + v — r_ 



(116) 



(117) 



where , r A , r + , r_ are the four real roots of A r 



For instance, for 
1 



ro,rs < r 



A 



one of the integrals we need to calculate is: 



A l dr 



(118) 

Indeed, we compute in closed-form: 



J a 



A l &7 



(r - r\)yj{r - a)(r - (3){r - yj(r - S) 



Pi 

,1113 r+-2q g - 7 r+ - 2a (3 - S r+ - 2a r+ - (3 r+ - 2a \ T(l/2) 
".2' '2' 2' '2'r+-2/?'a- 7 r+-2/3'a- ( 5r+-2/3'r+-ar+-2/3y T(3/2) 

(119) 



where 



Pi ■= 



- a r+ - 2/3 



(120) 



Also the radial integral involved on the LHS in the "balance" lens equation 
Q is computed exactly in terms of the hypcrgeometric function of Appcll F\ : 



rt/2 



dr 



\/|(CA + 3S 2 (1 + f ( a - $)2) L v/(r-a)(r-/3)(r-7)(r-(J) 



Pi 



r ( i /2) Fi (i^i^, 



\/£ V w (7"«)(^-«) r ( 3 / 2 ) V2'2'2 , 2 , a- 7 r+-2/3'a-5 



7 r| - 2a /? - <5 r| - 2a 
-2/? 



(121) 
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where £ := ±(QA + 3S 2 (1 + A( a _ $) 2 ), w := ^— | and a, 0,7,5 denote the 

roots of the quartic polynomial i? in the presence of A eqn©. 

A complete phenomcnological analysis of our exact solutions in the presence 
of the cosmological constant A will be a subject of a separate publication [15] . 
Nevertheless, it is evident from the closed form solutions we derived in this work 
that the cosmological constant does contributes to the gravitational bending of 
light. 

10 Conclusions 

In this work the precise analytic treatment of Kerr and Kerr-de Sitter black holes 
as gravitational lenses has been achieved. A full analytic strong-field calculation 
of the source, image positions and the resulting magnification factors has been 
performed. A full blend of important functions from Mathematical Analysis 
such as the Weierstrafi elliptic function p and the generalized multivariable 
hypergeometric functions of Appell-Lauricella Fjj were deployed in deriving 
the closed-form solution of the gravitational lens equations. From the exact 
solution of the radial and angular Abelian integrals which are involved in the 
lens equations we concluded the A does contribute to the gravitational bending 
of light. A full quantitative phenomcnological analysis of gravitational lensing 
by a Kerr deflector in the presence of A is beyond the scope of this work and 
will appear elsewhere [T3] . We provided examples of image-source configurations 
that solve the gravitational Kerr lens equations and exhibited their appearance 
on the observer's image plane as they will be detected by an equatorial observer 
(6>o = 7r/2,0 o = 0) and an observer located at 8o = 7r/3,</> = 0,for various 
values of the Kerr parameter a, and the first integrals of motion <£>, Q. 

The theory produced in this work based on the exact solution of the null 
geodesic equations of motion in Kerr spacetime wll have an important applica- 
tion to the Sgr A* galactic centre supermassivc black hole [15] . It may serve the 
important goal of probing general relativity at the strong field regime through 
the phenomenon of gravitational bending of light induced by the spacetime cur- 
vature. It is complementary to other investigations which have the ambition to 
probe gravitation at the strong-field regime through the relativistic effects of 
periastron precession and frame-dragging [T9] . 

There is a fruitful synergy of various fields of Science: general relativity, 
astronomy, cosmology and pure mathematics. 
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A Transformation properties of Lauricella's hy- 
pergeometric function Fp. 

In this appendix we prove useful transformation properties of Lauricella's hy- 
pcrgcometric function Fr>. We first introduce the function and its integral rep- 
resentation: 



Lauricella's 4 hypergeometric function of m- variables. 



Fd(q;,/3,7,z) 



E 



m ,ri2 — 



(oQni + -n m (/3i)n 1 '■■(j) 
(7)n 1 + ...+n m (l)n 1 " " (1) 



- 1 



(122) 



where 

The Pochhammcr symbol 
r(a + m) 



Z — (zi , . . . , Zra) , 

= (/?!,..., /3 m ). 



{oi)m = (a, m) 



is defined by 



(a)m = 



r(a) 



1, if m = 

a(a + 1) • • • (a + m - 1) if m = 1,2,3 



The series admits the following integral representation: 



(123) 



(124) 



Fu(a,/3,7,z) 



r( 7 ) 



r(a)r(7 - a) Jo 



t a -\i - ty-"- 1 ^ - Zl ty^ •••(!- z m t)-' 3 -dt 



which is valid for 



(125) 

Re(a) > 0, Re(7 — a) > 0. . It converges absolutely inside 

(126) 



the m-dimensional cuboid: 

N < 1,0' = 1, ■■■,m). 
Proposition 2 The following holds: 
F D (a,0,0',0",y,x,y,z) - (1 - z)- a x 



z — x z — y z 
z — 1 z — 1 z — 1 

(127) 
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Proof. Applying the tranformation in equation 

v v 



1 — z + vz 



(1-z) 



1 VZ 

1 z-1 



(128) 



we get: 



z-1 



1 - V- 



VZ 

z-1 



-P' 



, {t-uz)- 



VZ 1 
2-1 J 



(129) 



Thus 

IRF D = (1 - z y a x 



^(l - v y~ a -\\ - ,/£_£)-/> (i - v z -JL)-p\x _ I/ _!L)-(7-/»-/»'-/»") 

(130) 



and proposition follows. 



Proposition 3 The following identity holds: 



— 2p?L 2 Zm ~ z z 1 Fc (I i I _i 5 gdl I 23 ) z j ~ z s) Zj \ 

|o| y z m \J Zj - z 3 ^z m - z 3 \2' '2' 2' 2' Zj - z 3 ' z m (zj - z 3 ) ' Zj - z 3 J 

J_ 3 j~z£ J zm-Zs 1 x ( -£g F / 1 1 1 1 3 Zj(l - Z 3 ) Zj (z m - Z 3 ) Zj \ + 

\a\ V z ™ y 2j - ^3 y/zm - z 3 | 1 - z 3 \2 ' ' 2 ' 2 ' 2 ' zj - z 3 ' z m (zj - z 3 ) ' Zj - z 3 J 

1 F A 1 1 3 z 3 (z m - z 3 ) Zj \ 1 
1 - z 3 \ 2 ' 2 ' 2 ' 2 ' z m (zj - z 3 ) ' Zj - z 3 J J 

(131) 



Proof. We start with the integral representation of Lauricella"s hypcrgcometric 
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function: 



P f 1 1 1 13 



du u- 1/2 {l - uxi)~ l (l - ux 2 )~ 1/2 (l - ux 3 
r(3/2) f 1 du 1 1 1 - ux 3 



r(l/2) 7 ^/ul-Wi^T 1 !ii 2 \/l - ux 3 

r(3/2) 
r(i/2) 



(/(/ 



o v^(l - WXi) y/\ - ux 2 y/l - ux 3 



I 



du u x 3 



o \/u(l - uxi) $1 - ux 2 \J\ - UX~3 



w : + u „.„ — z 3 (l-z 3 ) — Zj (z m -z 3 ) 



B Time-delay assuming vanishing A. 



For the time-delay, in the case of vansihing cosmological constant, we derive the 
equation: 

rr r 2 (r 2 + a 2 ) , T 2GMr , 2 . , f a 2 cos 2 0d6 , s 

(133) 



ct = / — =^dr + / = (a -$a)dr + / 



±\/e 



In calculating the last angular term in (|133|) and using the variable z = cos 2 9, 
one of the integrals we need to calculate is: 

1 a 2 

(134) 



2 l a l Jo \fz(z m ~ z)(z- z 3 ) 
Indeed, its calculation in closed analytic form gave us the result 
1 a 2 [*i dz z 



2 \a\ J ^Jz{z m - z)(z - z 3 ) 

la 2 z 2 (z m - Z] ) p / 3 1 5 Zj_ z^z^_-z 3 \ T(l)r(3/2) 



2 M Zm^/fz, - z m )(z 3 - Zj)z~j V 2 2 2 z m z m Zj - z 3 J T(5/2) 

(135) 
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In total we derive for the angular integrals in (|133|) : 
r ° a 2 cos 2 6d0 



^time— delay 



la 2 (z m -z s )z m (1 113 z m -zs z m -z s \T(l/2) 



2H y/z m (z m -z s )(z m -z 3 ) V 2 ' 2'2'2' z m ' z m - z 3 7 T(3/2) 

h • /a M lfl2 z 2 s {z m -z s ) 

[1 - sign(6»s o e mS )\- 



2 M Zm\/( Z S - Z m )(z 3 ~ Z S )z S 



Z 1 5 z s z s z m -z 3 \ T(l)r(3/2) 

Xi< l i ' n> oi o! 1 ■ 



Z Z Z Z m Z m zs~z 3 J r(5/2) 
la 2 (z m -2o)^m c, / 1 1 1 3 z m - z z m -z fl \r(l/2) 



2|a|Vz m (z m -zo)(z m -z 3 ) \ 2 ' 2 ' 2 ' 2 ' ^ ' z m - z 3 J T(3/2) 
[1 - sign(0 o ° ^mo)]^ 



xFi 1, 



2 l a l 2 mV / (z - z m )(z 3 - z )z 
3 1 5 z z z m - z 3 \ r(l)r(3/2) 



Z Z Z 2 m Z m Zq 23 ^ / r(5/2) 

2 1«| Vz2 n (z m -z 3 ) r(2) V 2 2 '*»-3s/ 

We now turn our attention to the calculation of the radial contribution to 
time-delay in equation f|133[) . Indeed, the first term can be written: 

rs r 2 (r 2 +fl 2) 

dr 



/ 1 a 2 -2GMr\ 

r 2 dr f s 2GMr f rs 2a 2 GMrdr 4G 2 M 2 f rs I 1 A J 

— — + / dr - / 1 / -dr 

R Ja c 2 y/R J a c 2 Ay/R C 4 J a y/R 

(137) 
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In total for this radial term the exact integration yields the result: 



rs r 2 (r 2 +a 2 ), a 2 ft"z s 
dr = — r^=^= 

Ay/R j r s -a 



i 1 oo 1 1 3 f3r s -a rs- a r s - a /3 - 7 S - j3 rs - a\ T(l/2) 
xI, d -,-2,2, 



2' ' ' 2' 2 7 2' ar s - P' r s ~ r s - /3 a--/' 5 -ar s - (3 J T(3/2) 
ail"zs 2GM 1 1 1 1 3 /3 rs — a rs — a rs — a j3 — 7 1 rs — a 

^-D 77) _1 > !j 



r s ~a c 2 \ 2 '2 2'2'ars — /3'rs — /3'rs — /3a — 7'wrs — /3 

rs— 

z s Q"A* + d 



(r+ -a)yj 



.1113 (3-r+rs-a r s - a /3 - 7 <J - /3 r 5 - a^ T(l/2) 

x b D 77) X ) 77' 77' 77) ' 



2' ' 2' 2 ! 2' a-r+ r s r s -^a-7' S-ar s -l3j r(3/2) 

r s -a / 3 1 1 5 /? - r + r 5 - a r s - a /3 - 7 (5-/3 r s - a ^ T(3/2) 
J-^C 77> V 



r s -/3 V 2 ' ' 2 ' 2' 2' a-r+ r s - /3 ' r s -/3a- 7' <S - a r s - /3 / T(5/2) 



(r_ - a)y^ 



rg-a 
rs— ft 



1113 /3-r_r s -a r s - a /3 - 7 S - /3 r s - a\ T(l/2) 

x 77> ^ 77) 77) 77) ' 



2' ' 2' 2' 2' a - r_ r s -/3' r s - /3 a -7' 5- a rs - (3 J T(3/2) 
rg — a /3 1 1 5 /? — r_ rs — a rs — a /3 — 7 (5 — /3 ?*s — a\ T(3/2) 

T-f'D 77) 1) 



r s -/3 V 2 ' ' 2' 2' 2' a-r_ rs -,3' rs -/3a-7' <5- ar s -/3y T(5/2) 

4G 2 M 2 0"zs /l 1 1 3 rs-a^-7 5 - /3 r s - a\ T(l/2) 

? 1 l2'2'2'2'rs-/3a-7'(5-ars-/3/ T(3/2) 

V r s-0 

4GHP Atd il"(-zs) 



„4 "+1 



'1 1 1 1 3 /3 - r+ r s - a r 5 - a /? - 7 J - j3 r s - a\ T(l/2) 



2 ' ' 2 ' 2 ' 2 ' a - r+ r s - /3 ' r s - (3 a - 7 ' 8 - a r s - (3 J T(3/2) 
r s -a / 3 1 1 5 /3 - r + r s - a r s - a - 7 5 - /3 r s - a\ T(3/2) 

J-^B 77) X > 77) 77 ) 77) ' 



r s -/3 V 2 ' ' 2' 2' 2' a-r + r s -/3' r s -/3a- 7 ' 5-ars -^7 T(5/2) 



, I 'I 1 1 1 3 ^~ r - r ^" 



2' '2'2'2'a-r_rs-,3'rs-/3a-7'(5-ars-/?/ T(3/2) 

?'s — a /3 1 1 5 j3 — r_ rs — a rs — a /3 — 7 (5 — /3 ?*s — a\ T(3/2) 
7^ 77)1)77)77) 77)' 



r s -/3 V 2 ' ' 2' 2' 2' a-r_ r s -/3' rs-/3a- 7 ' 5-ars -/3/ T(5/2) 
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